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ABSTRACT 

We consider a generalized model of seismic-wave propagation that takes into account 
the effect of a central magnetic field in the Sun. We determine the g-mode spectrum in 
the perturbative magnetic field limit using a one-dimensional Magneto- Hydrodynamics 
(MHD) picture. We show that central magnetic fields of about 600-800 kG can displace 
the pure g-mode frequencies by about 1%, as hinted by the helioseismic interpretation 
of GOLF observations. 
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rS I 1 INTRODUCTION 



Currently there is very little direct information about the 
structure and strength of magnetic fields in the radia- 
tive zone (RZ) of the Sun, for a short rev iew see In- 
troduction of the paper iBurgess et al.l l2004al) . Some au- 
thors argue that for the young S un (< 3 — 30 Myr) 
relatively smaU fields, < 10 kG fMoss 2003) a nd < 
1 Gauss jKitchatinov. Jardine fc Collier-CameronI 1200 if) 
could survive, being relic fields captured from the primor- 
dial ones in the protostar plasma. For the Sun at the present 
epoch there is an upper bound of 2 — 3 MG near the 
tachocline obt ained from the magnetic sp litting of acous- 
tic oscillations jRuzmaikin fc Lindsevl2002ll . However, some 
authors have consid ered very strong magnetic fields in the 
RZ, up to 30 MG iCouvidat. Turck-Chieze fc KosovichevI 

liool). 

Here we suggest a new way to estimate the magnetic 
field strengths in the RZ of the Sun by relating them to 
the frequency shifts of g-mode candidates suggested by the 
first observations made with the GOLF (Global Oscillation s 
at Low Frequencies) experiment iXurck-Chieze et alJl20oJ) . 
We discuss some effects of RZ magnetic fields which could 
explain the displacement of g-mode frequencies with respect 
to the theoretical frequencies calculated in the absence of 
magnetic field. Indeed the existence of such shifts are hinted 
in golf's data. If eventually confirmed by further data, 
the idea that RZ magnetic fields cause such frequency shifts 
would provide us with a useful tool to estimate their mag- 
nitude. 
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In order to find spectra of seismic waves accounting for 
the magnetic field in the RZ a number of assumptions is 
required. For example: 

(i) We consider ideal MHD neglecting both the heat con- 
ductivity and viscosity contributions to energy losses, as well 
as the ohmic dissipation. 

(ii) We linearize the MHD equations about a static back- 
ground configuration, i. e. a background configuration which 
is time independent and for which the background fluid ve- 
locity vanishes, vq — 0. 

(iii) We assume the fluctuations to be adiabatic, with the 
contributions of fluctuations to the heat source vanishing: 
Q' = 0. 

(iv) Moreover, we consider a fully ionized ideal gas, so 
that the thermodynamic quantity, first adiabatic exponent 
7 = Cp/cv, is time independent and uniform. For numerical 
estimates we will take 7 = 5/3 for hydrogen plasma. 

(v) We adopt the Cowling approximation, which amounts 
to the neglect of perturbations of the gravitational potential, 
(i.e.: <t>' = 0). 

(vi) We assume a rectangular geometry with Cartesian 
coordinates: x, y, and z, where z corresponds to the so- 
lar radial direction. The background quantities vary along 
the z direction only (which implies the local gravitational 
acceleration, g, is directed along the z axis, but in oppo- 
site direction). We also take a constant, uniform background 
magnetic field. Bo, pointing along the x axis. 

(vii) The background mass-density profile is assumed to 
be exponential, po ~ pc exp[— z/i/], for constant pc and 
H. The conditions of hydrostatic equilibrium for the back- 
ground then determine the profiles of thermodynamic quan- 
tities, and in particular imply 7 is a constant. We assume 
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that the Brunt- Vaisala frequency is zero in the convective 
zone (CZ) and non-zero, but constant in the RZ. 

In what follows, we shall again specify the assumptions 
used, as they are needed, in order to keep clear which results 
rely on which assumptions. 

Note that, deep within the radiative zone, the last ap- 
proximation above holds to very good approximation for real 
mass-density profiles obtained by Standard Solar Models, 
provided we identify the 2 direction with the radial direc- 
tion. The constancy of 7 in this region is also expected since 
the highly-ionized plasma satisfies an ideal gas equation- 
of-state to good approximation. The rectangular geometry 
provides a reasonable approximation so long as we do not 
examine too close to the solar centre. What is important 
about our choice for Bo is that it is slowly varying in the 
region of interest, and it is perpendicular to both g and all 
background gradients, Vpp, Vpo, etc. ^^^^ 

As suggested in (jBurgess et al.l |2QQ4j such one- 
Dimensional (1-D) picture can be fully described in ana- 
lytical terms in contrast to the 3-D case. There are two 
param eters which describe the spectra of magneto-gravity 
waves iBurgess et al. 2004a') : (i) strength of the background 
magnetic field Bo and (ii) the dimensionless transversal wave 
number K — kxH . Here H is the density scale height and 
kx is the projection of the wavevector onto the x axis. Let 
us estimate the value of the transversal wave number that 
could be relevant for the g-mode candidates observed on the 
photosphere. 



Since g-modes decay in the CZ as 



-Kz/H 



only 



modes with low transversal wave number if ~ 1 — 4 (long 
wave lengths) could be seen at the photosphere. This follows 
from the simple estimate for the longitudinal fluid velocity 
Vz{z) which is directed along the Sun-Earth line and causes 
the Doppler shifts of optic lines registered by the GOLF 
experiment: 



NH b, 
K B, 



-3K o -1 

■e ~ 2 mm sec 







(1) 



This formula comes from Eq. (14) (equivalent to our 
Eq. Cnj) and Eq. (30) of Ref. l|Burgess et al.l l2004a^ for 
the decaying solution Bi^\z) = fe^e"^^^^, where B^^^ is 
the 2-component of the magnetic field perturbation. 

Here in the right hand side we substituted the sensi- 
tivity of the GOLF-instrument to the minimum fluid speed, 
Vz — 2 mm/sec, while in the left hand side we substituted 
the frequency estimate lj N and the wave length through 
the CZ: Rq — zrz ~ 3H — 0.3Rq. For instance, substi- 
tuting for the magnetic field perturbation, bz/Bo = 0.01, 
N = 2.8x 10"^ rad sec for the Brunt - Vaisala frequ ency 
in the RZ, H = O.IRq = 7 x 10^ cm llBahcalllllQsS) for 
the density scale height, one obtains e~^'^/K ~ 10"^ from 
which the estimate K — Kmax ~ 4 comes. 

We organize our presentation as follows. In Section |21 
we formulate the MHD model for an ideal plasma. In Sec- 
tion |3 we linearize the full set of MHD equations and then 
derive a single master equation for the z-component of the 
perturbation velocity v^^\z). This component leads to the 
Doppler shift of the optic frequencies measured in helioseis- 
mic experiments. In Subsection 13 . 1 1 we check the validity of 
the master equation against the well-known case of stan- 
dard helioseismology in an isotropic plasma, without mag- 



netic fields. In Subsection 13 . 21 we derive the simplified mas- 
ter equation in the perturbative limit. In this limit there 
are no MHD (slow or Alfven) r esonances w ithin the Sun, a 
situation which was treated in liBurgess et al. 2004a') . This 
perturbative method allows us to use a standard quantum 
mechanical 1-D approach to determine an exact analytical 
spectrum of g-modes in the presence of RZ magnetic fields. 
In Section^ we summarize our results. 



2 BASIC IDEAL MHD EQUATIONS 

We describe the Sun as an ideal hydrodynamical system 
characterized by nonlinear MHD equations. The mass con- 
servation law for the total density p can be written as 



|+V.(pv)=0. 



(2) 



The total pressure p is a sum of equilibrium, po, and non- 
equilibrium, p, parts, p = Po -I- p. Since viscosity can be 
neglected, momentum is conserved according to 

9v 



dt 



+ (v ■ V)v = --VP + -!-(B ■ V)B + g. 

p 47rp 



(3) 



Here the total pressure P = Po{z) + P consists of two terms: 
the equilibrium and non-equilibrium parts. The first is ex- 
pressed as Poiz) ~ Po{z) + Bq/Stt and obeys the equation 
VPo = po(z)g. The non-equilibrium part of the total pres- 
sure P =p + (l/47r)(Bo ■ B') + B' /Stt involves non-linear 
terms coming from the total magnetic field B = Bo + B'. 

The evolution of the magnetic field is governed by Fara- 
day's equation, 

^ = (B ■ V)v - (v V)B - Bt., (4) 

where u = V • v 7^ is the compressibility of the gas. Finally 
the conservation of entropy leads to the energy conservation 
law for an ideal plasma. 



d_ 

dt 



where p — po + p is the total gas pressure. 



(5) 



3 LINEAR MHD MASTER EQUATION 

For definiteness here we consider the same generalized he- 
lioseismic model already proposed in Ref. iBurgess et al.l 
l2004a|) , adopting an approximately rectangular rather than 
cylindrical, geometry. In this case it is convenient to choose 
a Cartesian coordinate system whose z-axis is the "ra- 
dial" direction; opposite to the local acceleration, g — 
(0, 0, —g{z)). With this choice we take z in the range (0, Rq), 
where z — represents the solar center and z = Rq de- 
notes the solar surface. The radiative zone corresponds to 
z < 0.7Rq. 

The model assumes the background magnetic field to be 
directed along the x-axis Bo = (-Bo(z), 0, 0), ensuring that 
physical gradients lie along the z-axis. Such a field mim- 
ics both a toroidal field lying in the equatorial plane and a 
poloidal field perpendicular to the equatorial plane. 

We now linearize the above MHD equations I2I5II . so 
that all variables are split into background and fluctuating 
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Euler quantities, f ~ fo + Z'^', with /^^•' denoting small 
fluctuations about the background value fo- In what fol- 
lows we also neglect rotation of the Sun, so that the back- 
ground velocity is zero, vo = 0. Moreover, we assume that 
the MHD perturbations are independent of y, which implies 



that B 



(1) 



0, leaving just six perturbation functions. 



instead of eight. 

I n addition, following ifRuderma n. HoUwee fc Goossenj 
Il997l) . we consider plane waves propagating along the x- 
axis and therefore the dependence of all functions on the 
coordinates {x, z, i) can be reduced to just 9 = x — Vt, with 
V = ui/kx being the phase velocity. This can be done since 
all functions depend only on the harmonic factor g»'=a:a;-ici't 

This way all perturbations are functions of two vari- 
ables, f^^\9,z), and the linearized MHD equations involve 
partial derivatives over 6 and z. Notice that d/dt — ~Vd/d6 
and d/dx — d/dO. 

Let us rewrite the initial system for the six functions 
f^^\B,z) in Eqs. (HEJ expressing all of them through just 
two quantities, v'"^\z,0) and P''^\z,9). One obtains: 



9P(1) _ poDa dv 



(1) 



dz 
dz 



V 09 

V ap(i) 



9Si 



(1) 



F i 
Dc 39 



(1) 



^'i^^ dpo pagDA (i) 
dz DcV • 

,.(1) 



89 



poDc 89 V dz ' 



,.(1) 



-Bo 



V ' 



Dc ' 



(6) 

(7) 
(8) 
(9) 
(10) 
(11) 



89 poDc 89 

We introduced above the important coefficients F, Da, Dc 
which in turn are functions of "z" via po{z), the Alfven ve- 
locity va{z), and the sound velocity Cs{z): 

PoDc 



Dc 



/2 - ci ' 



Da = V 



2 

- "A, 



4) 



(12) 

where = 'v'a(?sI{i'\ +Cs) is the squared cusp velocity. The 
zeros of the last coefficient Dc correspond to the slow [V = 
ct) or Alfven {V « va, va ^ c^) resonances respectively. 

Differentiating Eq. © over 9 and using the next two 
equations one gets the master equation for the function 



>{9,z) as: 



8z 8z ^ 



.2„,(1) 



PoDa 



8'v. 



89'^ 
d 



2 F dg 

9 TT+POT'+aT- 
Vc dz dz 



Fcj 
Dc 



„(i) 



(13) 



It is easy to see that, neglecting gravity, our generalized 
eq uations l|S)l-l[T^ recover the total system (34)- (38) derived 
in llRuderman. Hollweg fc Go ossenj| l997^ and Eq. I|13[l co- 
incid es with Eq. (39) i'n~ i^idBrmM^^ollweg fc Goossenj 
119971) . In case of the constant gravity, g=const, while 
the sound velocity Cg{z) = ^gH{z) depends on the 



var ying height scale H(z), E q. 1131 recovers Eq. (7) 
in jMiles. Allen fc Robertslfl992l) . 

For simplicity we consider below the particular case 
of uniform and constant magnetic field, _Bo=const, con- 
stant gravity, (7=const, and density scale height, _ff=const. 
Therefore the sound speed Cs and the Brunt- Vaisala fre- 



quency are also constants in RZ, 



7(;ff =const. 



7V^ = {g/H){'y — l)/7=const. This is somewhat of a sim- 
plification of the real Sun, but it allows us to derive qualita- 
tive estimates of the magnetic field corrections to the pure 
g-mode spectrum. 

In what follows in our numerical estimates we will use 
the following values: 7 = 5/3 (hydrogen plas ma), N = 2.8 x 
10~^ rad sec-\ H = O.IRq = 7 x 10^ cm (lBahcalllll98ll . 
hence g = 1.37 x 10^ cm s~^, Cs = 4 x lO'^ cm/s. 

Separating in master equation 11311 the exponential de- 
pendence e''^^* and taking into account the background den- 
sity profile po = pce~^^^ , one gets the ordinary second order 
equation^ : 



1-^11 



2 (1) 



d^v 



(1) 



dz2 



1 dv 
H~dz 



+ 



V2 



„(1) 



: 0. 



(14) 



This equation generalizes Eq. (16) of our 
work (Burgess et al. 2004a) accounting for the compress- 
ibility parameter ai = V^/cl = uj'^/klcl{= 0.24uj'^ / K'^ N'^ 
for hydrogen plasma with 7 — 5/3). Note that ai is small 
for very low frequencies of g-modes cu <ti N, ai <^ 1, 
and therefore it was neglected in the problem considered 
in (iBurgess et al...2004a'l . 



3.1 Zero magnetic field limit 

For VA = Bo — with the use of the transformation 
v'-J-^z) = *(z)e^''^-^ Eq. lO converts into 



di2 



+ 



4^2 



* = 0. 



(15) 



This 1-D MH D equation coincides with th e 3-D oscillation 
Eq. (7.90) in dChristensen-DalsgaardlboOSh : 



d^* 

dr2 



+ 



4//2 



+ kl 



N^{r) 



- 1 



* = 0, 



(16) 



here we have substituted the acoustical cut-off frequency 
LOc = (c?/4H2)(1 - 2dH/dr) ct/A H^ and put Sf =1(1 + 
l)Cs/''^ — fcfeCa, see Eq. (4.60) in lIChristensen-DalsgaardI 
2003). 

We would like to stress that, in the JWKB approxima- 
tion for low frequency g-modes, ^ V^, Sf ^ lj^, both 
equations 1151 and II6I I lead to the analogous spectra: 



U!g{n, kx) = 



7r(n-hl/2) 



ZRz (ID model). 



u)g{n,l) 



./lilTT) j;:^- N{r)dr/r 
Tv{n + l/2 + ag) 



(3D model), (17) 



^ Here we have corrected the factor in front of the first derivative 
in Eq. (16) of jSurgcss et al. 2004a) changing (7 — l)/7_ff — > 1/H. 
This does not affect the MHD spectra found there. 
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Figure 1. The g-mode frequencies for n = —1, —2, . . . , —10 nor- 
malized to Brunt- Vaisala frequency, co/N, versus wave number, 
K. Solid lines are for zero magnetic field, while dashed lines cor- 
respond to a fixed magnetic field Bq = 700kG. 
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Figure 2. The resonance positions are shown in the plane Bq — K 
(magnetic field vs wave number along magnetic field) . The curves 
separate the regions where perturbative or non-perturbative ap- 
proaches may or should be used, (^z = 1 means that at place z 
the Alfven resonance occurs. 



where n is the radial order of the wave, and the transversal 
wave number is the 1-D analogue of the angular degree I 
in the 3-D case. 

Let us denote the factor in brackets in Eq. 11511 as 



0- 



jV2 

7J 



(18) 



here K = k^H is the dimensionless wave number, K ^ 
Kmin ^ 2-kH/Rq ~ 1. Eq. H15|l can be rewritten using 
Eq. (|THJ as 



d;7 + ife*W = o- 



(19) 



In RZ where Nrz = N — const 7^ one obtains the solution 
of Eq. 1^ in the form *flz(z) = Crz sm{l3z/2H), where 
Crz is a constant. It accounts for the boundary condition 
Vz^\o) = at the center of the Sun. In CZ we approximate 
Brunt- Vaisala frequency by the value Ncz = 0. In that case 
(3 tends to iV — iy^4:K^{l — ai) + 1. Taking into account 
a second boundary condition, that there are no solutions 
which grows with z, one gets the decaying MHD wave in 
the form: ^cziz) = Ccz exp{-zr/2H). 

Now matching both solutions at the top of RZ, 
z — zrz, (for log a rithm ic derivatives see the book 
by iLandau fc Lifshit j ll2000D l. one obtains the dispersion 
equation in the case Bo = 0: 



/3cot 



PZRZ 



2H 



= -r. 



(20) 



or the g-mode spectrum in our 1-D model is given by 
(3zRz/2H + TTn = arccot[— r//3], n = —1,-2, —3,... (see 
solid curves in the Fig. 1). 



3.2 Magnetic corrections to the g-mode spectrum 

In order to obtain spectra of g-modes in the presence of 
magnetic field let us define the coefficient in front of the 
second order derivative in Eq. 1141 as 1 — where ( = 
'"Ao^'^^^i^ ~ ai)/V^. The Alfven velocity at the solar center 



is given by vao = Bq / A-k pc ~ (i3o/43.4 Gauss) cm sec~^, 
here, and through the paper, pc — 150 g/cm"^. 

We consider the perturbative regime for magnetic fields, 
where vl = v\oe'^" < = cj^kl = [N"" K^'jicu/Nf , 
so that C ^ 1 fl-nd MHD resonances {(, = 1) do not appear 
within the RZ. This region is the lower one in Fig. |21 The 
dashed curves labelled ^0.3 = 1 and (^0.7 = 1 correspond 
to resonances that occur at O.3-R0 and 0.7Rq in the non- 
perturbative region. The solid curve is chosen to illustrate 
the separation between the two regimes, according to the 
criterion ^0.7 ~ C(O.77?0) — 0.1 ^ 1. The perturbative mag- 
netic field for which the maximum Alfven velocity is small, 
v\{0.7Rq) < V^ obeys 

2.8 X lO'^ r. ,r ^ ATN 

Bo < Gauss (for a; < N). 

The absence of MHD resonances allows us to set 

-1/2 

Z 

e 



(21) 



,,(1) 



ii 

\/2 



r.2 



and thus derive from Eq. 11411 the 0(i;^/V^)-correction to 
Eq. Cnj: 



d^* 

d^2 



+ 



4H^ 



k% 



*(^) =0, 



(22) 



here the wave number correction kg is given by 

"iV2 



7 2 t^ ,2 

Ks = (1 - ai)y^k^ 



+ ai 



Introducing the notation: 



>0, 



(23) 



one gets from Eq. 12211 : 



dz2 



+ 



4^2 



-I- k^Age 



■■/H 



*(2) = 



and after the change 2s 
takes the form 



z/H + ln(4iC^Ao) this equation 
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^ + [/3^ + e^1*(.) = 0. 



(24) 



The general solution of Eq. 1241 is expressed via the Bessel 
funct ions of the fi rst kind, '^(s) = CiJipie") + C2J-ii3{e'') 
(see iGradstein fc Rv zhik (2000)), Ci and C2 are constants. 
Now, coming back to the variable z, we use the boundary 
condition at the solar center, ^'(0) = {vz{0) = 0), to obtain 
the solution of Eq. II24II in RZ, ^ z ^ zrz, 



vE'hz(z) = Crz 
Ji^{2KAl^ ^] 



1/2 z/2Hx 



1/2 



1/2 z/2H^ 



(25) 



In CZ the Brunt- Vaisala frequency A'' vanishes, A'^ = 
0, so that we are led to the same solution as in isotropic 
case ^cz[z) = Ccz &qi{—zT / H), neglecting the possible 
existence of a CZ magnetic field (if such a field is present, 
with strength of 300 kG, the g-mode frequencies change by 
no more than 10~^). Since the arguments of Bessel functions 
in Eq. (1241 are small, we can use the first terms in the Bessel 
function series only. 



J. (2) 



Z. fc!r(. + fc + i) I2J ' l^'^f-K-' 



fc=0 

1/2 



J±,li{2KA^,''e 
l±il3 



{KA^ 



1/2 z/2i^^±^/3 



r(i±i/3) 



1 - 



+ 0{K^Aoe 



z/H. 



(26) 



Then by matching the logarithmic derivatives of the solu- 
tions 'i! RZ and "^cz at the top of RZ, z — zrz, one obtains 
the generalized dispersion equation for the case Bo 7^ 0: 



P[ 1 + 



X cot 



o 2 



2H 



1-13 



'2kr.z 



(l + (3^)sm{f3zRz/H) 



2Ki 



= -r, 



(27) 



z/2H 



where krz ~ KA^^^e^ 

This is our main equation. Clearly, when the magnetic 
field correction tends to zero, hence 2Kjiz/{^ + /3^) tends to 
0, one recovers the isotropic case, given in Eq. 1201 . 

We search for a perturbative solution of Eq. 1271 of the 
form P = (3o{1 + 5b) and w = ^Jo(^^)(l+aB(^^•)) where (3 = Po 
and correspondingly u — loq (n) are the solution of Eq. 12U1 
for Bo — and the smallness of 5s ^ 1 and as(ri) 1 
follows from ^fj^ -ti 1. 

The g-mode spectra for fixed Bo = 700 kC are shown on 
Fig. (see dashed curves). Figs.|3|and 2]display our results 
for as as a function of mode radial order n, magnetic field 
and wave number K. One sees that the magnetic field shift 
as of a g-mode frequency uJo{n) is always positive, qs > 0. 
In Fig.|3]we show the absolute values of the shift qs(w) for 
different g-modes and for fixed Bo = 700 kC as a function of 
K. Conversely, in Fig.^Jwe fix the wave number K = 2 and 
plot as(n) as a function of Bq. One sees that, the higher the 
mode radial order \n\, the less the magnetic field strength 
required to produce a given g-mode frequency shift. 



0.015 



0.01 - 



0.005 




Figure 3. The shift of g-mode frequency as versus wave num- 
ber, K, shown for fixed magnetic field strength, Bq = 700 kG and 
different mode radial orders, n = —1, —2, . . . , —10. 




1x10 2x10 3x10 4x10 5x10 6x10 7x10 8x10 

[Gauss] 

Figure 4. The shift of g-mode frequency ag versus magnetic 
field strength, Bq, shown for fixed wave number, K = 2, and 
different mode radial orders, n = —1, —2, . . . , —10. 



4 DISCUSSION 

We have given a generalization of helioseismology to ac- 
count for the presence of central magnetic fields in the Sun. 
We have determined the resulting g-mode spectrum within 
the framework of a perturbative one-dimensional Magneto- 
Hydrodynamics model. 

There are three factors influencing g-mode observation 
in helioseismic experiments. First, such g-modes should have 
long wave lengths to penetrate the CZ: in our case K ^ 4, 
and in the 3D-model, low values of I. Second, the radial order 
n should also be low, otherwise, such low frequency g-modes 
are much below the present experimental sensitivity. Third, 
there is a strong influence of the RZ magnetic field. 

If the magnetic field is too strong (more than a few 
MG ) all g-modes are lo cked within the Alfven cavity 
(see lBurgess et all (|22Q43))i hence these g-modes decay far 
beneath the CZ becoming invisible in the photosphere. This 
happens because the MHD energy in the radial direction is 
fully diverted to the transversal plane at the Alfven reso- 
nant layer position. On the other hand, for very small mag- 
netic fields (Bo ^ 1 MG) the magnetic field effect is neg- 
ligible and can not account for a possible discrepancy be- 
tween present experimental data and theoretical predictions 
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of seismic models. In contrast, we have shown that a solar 
radiative zone magnetic fields of intermediate magnitude, of 
the order 600-800 kG, can displace the pure g-mode frequen- 
cies by about 1% with respect to our model of seismic wave 
propagation, a value close to what is hinted by results of the 
GOLF experiment. 

We find that higher modes require a smaller magnetic 
field to produce a given g-mode frequency shift. However 
encouraging this result may sound, let us stress again that 
in our simple one-dimensional MHD picture (with further 
assumptions, such as adiabatic oscillations, exponential den- 
sity profile, constant gravity, etc) we can only make a quali- 
tative estimate of the magnetic field corrections to the pure 
g-mode spectrum. For example, our formalism can not ex- 
plain the magnetic splitting of g-mode frequencies over az- 
imuthal number m, as it requires 3-D. Further work in 3-D 
geometry is necessary to perform a quantitative comparison 
with t he frequency patterns ob served in the GOLF experi- 
ment jTurck-Chieze et al .120041) . Even within the simple an- 
alytic approach 1-D MHD model one may include viscosity 
effects in non-ideal plasma with finite conductivity, and also 
take into account magnetic field diffusion stabilizing MHD 
instabilities. 

Last, but not least, recall that our perturbative analy- 
sis avoids the appearance of MHD resonances that could 
lead to density spikes. These are potentially important, 
as they can affect neu t rino p ropaga tion through the so- 
lar RZ (iBureess et alJ |200A l2004allbl) . Improved deter- 
minatio n of neutrino mi xing parameters, e.g. by Kam- 
LAND iEguchi et al Eooi), allows one to carry out neutrino 
tomography in deep solar interior. Both regimes of "magne- 
tized helioseismology" (i) the MHD seismic models and (ii) 
the analysis of MSW neutrino oscillations in noisy Sun are 
complementary tools to explore RZ magnetic fields. A fully 
quantitative analysis may require the inclusion of th e differ- 
ential rotation in the RZ jTurck-Chieze et aljlioosll as well 
as non-linearities. 

Note added: As this paper was being typed, we saw 
a paper by Hasa n. Zahn. fc Christensen-DalsgaardI ll2005ll 
( astro-ph/05f f 472 ), where a similar idea to probe the in- 
ternal magnetic field of slowly pulsating B-stars through g 
modes is given. Their 3-D results are consistent with our 
simpler 1-D estimates. 
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